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QhI Introduction 

Oh: 

^ . Nonperturbative effects of non-abelian gauge theories take place in the infrared region. 

, Among such nonperturbative effects are the existence of the QCD vacuum with gluon and 

■ quark condensates, chiral symmetry breaking, confinement and the hadronization process. 
They can only be studied by approximate methods and in the framework of various effective 
models. For instance, one of the possibilities to approximately describe the gluon condensate 
is to introduce background color fields of certain configurations (see, e.g., 

Another example of nonperturbative problems is the physics of light mesons that can be 
described by effective four-fermion models such as the Nambu-Jona-Lasinio (NJL) quark 
model, which was successfully used to implement the ideas of dynamical chiral symmetry 
breaking {DxSB) and bosonization (see e.g. j^,j4| and references therein; for a review of 
(2+l)-dimensional four-quark effective models see 5[). 

In the framework of four-fermion models, it was shown that a constant magnetic field 
induces the dynamical chiral symmetry breaking (DxSB), as well as the fermion mass 
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generation, even under conditions when the interaction between fermions is weak. This 
phenomenon of magnetic catalysis was explained basing upon the idea of effective reduction 
of space dimensionality in the presence of a strong external magnetic field Q] (see also paper 
0] and references therein). It was also demonstrated that a strong chromomagnetic field 
catalyzes DxSB As was shown in jl|J], this effect can be understood in the framework 
of the dimensional reduction mechanism as well, and it does not depend on the particular 
form of the constant chromomagnetic field configuration. 

One of the solutions of the Yang-Mills equations that can serve as a model for the gluon 
condensate is a constant chromomagnetic field. Its role was demonstrated in jlll . 1^ . In these 
papers, the authors calculated the one-loop effective potential for a constant chromomagnetic 
field B and they demonstrated that it reaches its minimum at a nonvanishing value of B. 

However this simple analytical model of the gluon condensate with a uniform chromo- 
magnetic field B = const (the so-called "colour ferromagnetic state" ) ll|, |l2| suffered from 
an instability . This problem later has been studied in a number of papers. In particular, 
various methods were proposed to stabilize the solution by introducing a nonzero charged 
component of the gauge field. Various attempts have also been made to improve this by 
assuming a certain domain-like [lj| or non-abelian structures Il5| of the condensate field. 

nn 

Possibilities of thermalizing the system [16|, ll7[ and introdu cing a condensate of the time 
component of the gauge field Aq ^ were also considered |l[ llj 

The principal difficulty in finding a stable color ferrromagnetic state is that a local min- 
imum of the action can not be obtained, since the corresponding field configuration proved 
to be spatially inhomogeneous. In order to circumvent this difficulty, the method earlier 
applied in analyzing the quantum Hall effect [3] was used in 0, It was demonstrated 
that a spatially homogeneous state of the gluon field can be obtained by effectively reducing 



the dimensions of the prob lem from D = (3 + 1) to D = (2 + 1) and employing for gluons 
the technique used in [lij . for the description of the quantum Hall effect in a (2 + 1)- 
dimensional Fermi system. The resulting gauge configuration breaks the color symmetry 
and the arising color charge can be compensated by fermions (quarks) in the system inter- 
acting with the color field. It should be mentioned that in this case the quark density must 
be higher than in the hadronic phase of matter. 

At the same time, in the system with high fermion density, another nonperturbative 
phenomenon can take place, i.e., the color superconductivity (CSC) (see, e.g.|23^ - 
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Briefly, the mechanism of CSC can be explained as follows. Quarks in a dense medium 
become asymptotically free, their interaction is rather weak and attraction may arise between 
them. Components with longer interaction range are screened. As a result, quarks of 
different colors and flavors with opposite momenta may compose Cooper pairs, as electrons 
do in the superconducting metal, thus leading to the energy gain. One may expect that, 
similar to the case of the quark condensate, the process of diquark condensation can be 
catalyzed by intensive external (vacuum) gauge fields. For a (2+1)- dimensional model, 
this was recently discussed in jsQl. It is clear that a medium with the density ~ 10^^ g/cm^ 
needed for color superconductivity to arise cannot be produced in the laboratory. This dense 
media may be found, e.g., in compact stars [sil Is^. 

In the present publication, we investigate color ferromagnetism (CFM) and color super- 
conductivity as two mutually excluding possible phases in an SU{2) gauge field model. Our 
aim is to study in more detail the effective dimensional reduction in the gluon sector leading 
to stabilization of the ferromagnetic state, and also to consider for the role of fermions in 
stabilizing the system and their contributions to thermodynamical quantities in the color 
ferromagnetic phase. We analyse the ferromagnetic phase and compare the energy gains due 
to either of possible effects: creation of color ferromagnetism and color superconductivity. 

The paper is organized as follows. Section 1 contains the investigation of the gluon 
sector. Section 2 contains the discussion of the quark sector. In Sections 3 and 4, we 
consider the phase transitions at zero and finite temperatures respectively, and Section 5 
contains conclusions. 



I. THE EFFECTIVE DIMENSIONAL REDUCTION IN THE GLUON SECTOR 

Let us consider the gluon Lagrangian in the Yang-Mills theory with the SU{2)c color 
group for the field (a = 1,2,3) in the (3+l)-dimensional space-time 

A = -^(^;.)' - Y^iDtQ'.r+xaiD'uxb. (1) 

Here = + Q^, is the background field, are quantum fluctuations of the gluon 
field, = 5^^d^ — ge°'^'^X^^ is the covariant derivative in the background field, x? X ^'^^ 
ghost fields and {D^)ab = D'^f^D'^^, g denotes the gluon coupling constant (the Feynman 
gauge ^ = 1 will be adopted in what follows). Moreover, summation over repeated color 
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indices a — 1,2, 3; and Lorentz indices // = 0, 1, 2, 3 is implied. The above Lagrangian can 
be rewritten in the form demonstrating interaction of the colored ("charged") vector field 
with a — 1,2 and a "neutral", a = 3, i.e., "electromagnetic" field: 

^9 = -\fl - l\(D,W^ - D.W,)\' - igf,.W;W^ + ^^{W^W^ - W:^W,f, (2) 

with the following notations adopted: ^4^ = is the "neutral electromagnetic" field, 
fixv — d^Aj, — djyA^ is the corresponding field tensor; — ;^(^^ + iA^^) is the "charged" 
field, Dfj, = + igA^. We omitted the ghost contribution in the above lagrangian, having 
in mind that it will cancel the unphysical degrees of freedom in subsequent calculations. 

Let us assume that the background field is a constant abelian chromomagnetic field 
directed along the axis, playing the role of the "electromagnetic" field in the above 
equation. Choosing an asymmetric gauge for this field, we can write 

_ _ (3) 

In the linear approximation, the field equation for quantum fiuctuations of the color field 

in the background chromomagnetic field has the form: 

dIw^ - igTw, = 0. (4) 

We may write the following stationary solutions for two physical degrees of freedom of the 
fiuctuation field Wi^2 interacting with the chromomagnetic field B: 



W^,2 = e-'^^°+^'=^^^+^^^^^(ii„)i,2(xi - k^/gB), 



where the functions (wn) 1,2 (2^1 — gB) are expressed in terms of the known wave functions 
of the harmonic oscillator. The energy spectrum has the form: 

£2 = kl + 2gB{n + ^) + 2agB, (5) 

where n — 0,1,2,... is the Landau quantum number. It is clear that the energy might 
become imaginary, when the spin projection is a — —1, the Landau quantum number is 
n — Q and the longitudinal momentum projection is restricted by kl < gB (the so called 
tachyonic mode). This behavior of the solution demonstrates that the linear approximation 



for solving the equation for this particular mode is not valid Therefore, for this unstable 
tachyonic mode, which is denoted by the field 99, 



the exact nonlinear Lagrangian 



i^tach = \D^ip\^ + 2gB\ip\^ - ^\^\\ (7) 



2l|, 



should be considered. Unfortunately, the procedure of finding a nontrivial uniform vacuum 
state (/^ 7^ independent of x, as it was done in the case of the Higgs model, is not applicab 
here due to the presence of A^{x) in the covariant derivatives = df^ + igAfj^{x). In 
a new approach was proposed to achieve the stable vacuum state. There, the solution with 
^3 = 0, which is uniform in the X3 direction, was considered. This choice leads to stabilization 
of the solution, if the effective dimensional reduction D = (3 + 1) ^ D = (2 + 1) is 
implemented. 

It should be emphasized that the solution found in [2^ is stable only with respect to 
(2+l)-dimensional perturbations. However, the true vacuum state should be stable with 
respect to (3+l)-dimensional perturbations as well, i.e., against non-uniform perturbations 
along the X3 axis parallel to the chromomagnetic field direction. In fact, this is the condition 
that there should be no tachyonic modes with nonzero momenta along this axis, ^ 0. 

In order to find this condition, let us consider the unstable modes in more detail. It follows 
from (0) that unstable modes appear when n = 0, and the sign of the gluon spin projection 
is negative, while the longitudinal momentum is comparatively small, i.e., fcg < gB. In order 
to exclude tachyonic modes with ^3 7^ 0, the extention of the chromomagnetic field along the 
X3 axis L3 should be finite. In this case, inplying a periodicity condition for the "charged" 
field along this axis W^''^{x3 = 0) = W^''^{x3 = L3), the momentum becomes discrete: 

= 2Tm^/L^, where = 0, ±1, ±2, .... Thus, in order to exclude tachyonic modes with 
1^31 > 0, one should demand that the lowest non-uniform modes, with n-^ = ±1, should 
already have real energy, i.e., that they were not tachyonic, and this supplies a restriction 
on the maximum extension L3 of the chromomagnetic field B along the X3 axis, 

L3 < LT^ = (8) 

This limitation provides a physical meaning for the choice of the solution uniform along the 
0:3 axis, which implies a physically reasonable restriction, i.e., the field should exist only 
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inside certain domains with finite dimensions. In other words, when the length L3 (the 0:3 
dimension of the domain) is fixed, the inequahty (jH)) may be considered as a hmitation for 
the maximum value of the stable chromomagnetic field inside the domain 

9B < 9B0 = i^r. (9) 

Thus, we should not consider a solution, which is uniform in the whole space, but rather 
a single domain, inside which the field is constant in magnitude and in direction. This 
limitation for its maximum extension and the boundary conditions imply that the domain is 
surrounded by other domains with fields inside having different strength and orientation in 
configuration and group spaces, restoring the symmetries of the whole system. In fact, the 
above limitation for the L3 extension of the chromomagnetic field domain implies effective 
dimensional reduction for the unstable mode that results in its stabilization. As it was 
mentioned above, we restrict ourself by the simplest assumption of periodicity of the solutions 
along the X3 axis. The detailed investigation of the structure of domains and boundary effects 
that inevitably occure, when domains are put together, is out of the scope of the present 
article. 



II. QUARK SECTOR 

The gauge configuration that consists of a constant chromomagnetic field B and a uniform 
gluon condensate W^'"^ 7^ is color charged. Following [2^, one has to introduce fermions 
with finite density into the system that will interact with the gauge field configuration and 
preserve its color neutrality. In the case of larger quark matter density, an effective model 
describing the interaction of quarks should be employed. Let us first give some arguments for 
the choice of the QCD-motivated extended Nambu-Jona-Lasinio effective model of quarks 
introduced below. To this end, consider two-flavor QCD with nonzero chemical potential 
fi and color group SUc{2). In the previous section, we decomposed the gluon field Al{x) 
into a condensate background ("external") field A^{x) and the quantum fluctuations Q'^{x) 
around it, i.e. A^ix) = A"(x) + Ql{x). By integrating in the generating functional of QCD 
over quantum fluctuations Q'^{x) and further "approximating" the nonperturbative gluon 
propagator by a 5— function, one arrives at an effective local chiral four-quark interaction of 
the NJL type describing low energy hadron physics in the presence of a gluon condensate. 
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Finally, by performing a Fierz transformation of the interaction term, one obtains a four- 
fermionic model with (gg)-and (gg)-interactions and an external condensate field A^^{x) of 
the color group SUc{2) given by the following Lagrangian 

A = liYi^d, + gX{x)^) + /i7°]g + ^limY + (g«7'rg)2] + 
G 

+ Jf [iqc£e-f^q] [iqee-f^Qc] ■ (10) 

In (fTUj) . qc = C(f, (jc = q*C are charge-conjugated spinors, and C = ^7^7° is the charge 
conjugation matrix {t denotes the transposition operation), r = (r^,r^,r^) are Pauli ma- 
trices in the flavor space; {eY^ = e^^, (e)°^ = e"'^ are totally antisymmetric tensors in the 
flavor and color spaces, respectively. Clearly, the Lagrangian (ITUj) is invariant under the 
chiral SU{2)l x SU{2)r and color 5f/,(2) groups. 

As it was mentioned above, the physical vacuum of QCD may be interpreted as a region 
splitted into an infinite number of domains with macroscopic extension. Inside each such 
domain, there can be excited a homogeneous background chromomagnetic field p^], which 
generates a nonzero gluon condensate {FF) 7^ 0. Averaging over all domains results in 
restoration of the color as well as Lorentz symmetries 36|. 

Upon performing the usual bosonization procedure and introducing meson and 

diquark fields a, vr and 5, 5*, the four-quark terms are replaced by Yukawa interactions of 
quarks with these fields, and the Lagrangian takes the following form (our notations refer 
to four-dimensional Euclidean space with it = X4) 

A = -qihu^u + ifilo + a + ■i-f^Tn)q " ^(^^ + " ~ 

- S*[iq'Cee-f^q] - 6[iqee^^Cq% (11) 

where = d^ — igX^^aa/2 is the covariant derivative of quark fields in the background field. 

In order to investigate possible phase transitions in the quark matter in the framework 
of the initial model (jllll . we evaluate the path integral over meson and diquark fields by 
using the saddle point approximation, neglecting field fluctuations around the mean-field 
(classical) values < cr >= (Tq = 0, < vr >= ttq = and < 5 >= Aq, < S* >= Aq. 

Within this approximation, we obtain the quark contribution to the partition function 



Zq = exp We = / dqdq exp 



(12) 



8 



where 

A = -g(^7.V.)g - A*[zg*C£e37'g] - ^^[iqee^l^'Cq^ (13) 

with We being the Euchdean effective action, and Cq the quark Lagrangian. 

The partition function (fT^ is calculated in the standard way (for details see, js^). In 
principle, the gap Aq is complex. However, the partition function is real and depends only 
on the module squared of the gap. Its phase characterizes just the degeneracy of the vacuum 
and may be set here equal to zero. In this sense, it is understood that the following equations 
are expressed directly in terms of the module |Ao|, i.e. 

2Ao ^ |2Ao| = A. (14) 

In the present study, the background field is assumed constant and homogeneous, 
F"^^ = const. Then the Dirac equation 

(7^V^)^ = (15) 

for a quark with flavor i has stationary solutions ipk,i with the energy spectrum Ek^i, where k 
stands for the quantum numbers of the quark in the background field. In this case we arrive 
at the following Euclidean effective action: 

= ^ / ^ E (pI + ^" + - '^Z^)^) ■ (16) 

k,i,K 

Here, k = ±1 corresponds to charge conjugate contributions of quarks with color indices 
a = 1,2 (included in the quantum number k) and the spectrum j, moving in the 
background color field -F^^. Clearly, for a vanishing external field (-F" = 0), we have 



2 -Q 

4 = p ■ 



In the case of finite temperature T = 1/(3 > 0, the thermodynamic potential 
= -WeHPL^) [1] is obtained after substituting ^4 ^ ^(/ + i), / = 0, ±1, ±2, 



Nf l=+oo 
i=l K l=—oo k 



2n{l + 1/2) 



A^ + {ek,i - Kfif 



(17) 



Next, with the use of the proper time representation we obtain for the quark thermody- 
namic potential 



i—l K -^1^2 '—1 



As 



■)(-iy 



X 
k 



^exp[-s(A2 + (£;,,, -/s:^)2)], (18) 
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In order to find the phase of the system at zero and finite temperatures, solutions and 
energy spectrum e^^i of the Dirac equation (fT3jl in the chromomagnetic field should be used. 
For the case of a constant abelian chromomagnetic field, equation is decomposed into 
independent equations for quarks of different colors. Fermions are in the states determined 
by the following quantum numbers: color, spin, sign of the energy. Landau number and two 
components of the momentum p2, Ps- The corresponding energy spectrum is well known 

4,a,p, = 9H{n+^ + ^)+pl, (19) 

where o" = ±1 is the spin projection on the external field direction, ps is the longitudinal 
component of the quark momentum (— oo < p^ < oo), 

pl=gHin + ^) (20) 

is the transversal component squared of the quark momentum, and n = 0,1,2,... is the 
Landau quantum number. It should be mentioned that the boson condensate of the gluon 
field (f) may slightly change the structure of the fermion levels, and besides it may mix quarks 
of different colors. The former may be neglected as it will only slightly influence the total 
energy of fermions, and the latter mixing of colored fermions is just the mechanism that 
provides the ground state of the quark system with a nonvanishing color charge. 



III. CSC-CFM PHASE TRANSITION 



In electrodynamics. Cooper pairs can be produced only in the region, where there is no 
magnetic field, as the charged condensate forces the magnetic field outside the sample of 
a superconducting metal due to the Meissner effect. In the non-abelian case the situation 
is somewhat different. For instance, in the theory with the SU{2)c group, quarks in the 
fundamental representation form a doublet qi {i = 1,2), and the diquark condensate has 
the structure like < e^^qiqj >, which is a scalar in the color space. Hence, no contradiction 
between possible superconducting state creation and the presence of a chromomagnetic field 
arises 3J|. In the real QCD case with the SU{?))c group, the condensate has the form 
< ^^"^lalp >, where a,/3,7 = 1,2,3, and it is no more color neutral. This condensate 
expells the part of the color field, which might interact with it. In the group space, a neutral 
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chromomagnetic field can be written in the form of a superposition of commuting generators 
A3 and Ag, i.e., F12 = 51A3 + -B2A8. Therefore, for instance, the field F12 ~ A3 can in principle 
coexist only with the condensate of the form < e^^^^q^qp >. 

The SU{3)c group has a maximal abelian subgroup f/(l) x f/(l). In each of the subgroups 
f/(l) either color ferromagnetism or superconductivity is possible. In the present publication, 
we consider the SU{2)c gauge group. However, having in mind possible generalization for the 
SU (3)c group, we shall consider the situation, when color ferromagnetic and superconducting 
phases can exist only separately, and a phase transition can take place between them. 

First consider the case of zero temperature T = 0. In order to consider which of two 

phases, color ferromagnetic or supeconducting, is preferable, one should compare the energies 

of these two phases. The ferromagnetic phase is described by a constant chromomagnetic 
3 3 

field 

— ~^2i — const — B and the quark matter interacting with this field. The 
thermodynamic potential of this system is the sum of the chromomagnetic field energy 
5V(87r) and the total energy of quarks Eq that occupy the lowest Landau levels in this 
field: 

Eq = f2g|T=0,A=0 = ^qO, (21) 

where i^g|r=o,A=o is determined by the corresponding limit of formula ()18p (the first term 
in square brackets, and A = in the exponent). The contribution of gluon quantum 
fluctuations around the chromomagnetic field can be neglected assuming the value of the 
coupling constant is small, g'^/lAn) ^ 1, in the limit of larger quark density. 

In our model, the chromomagnetic field strength B and its extention L3 (the dimension 
of the domain along the chromomagnetic field direction) are related by the inequality (jH)) 
that determines the stability condition and defines the maximum possible chromomagnetic 
field strength Bq (jHl). We will make a numerical analysis of the thermodynamic potential 
of quarks and, hence, it is natural to introduce dimensionless variables a,x, and h{x,a) 
in place of the chromomagnetic field B, fermion density p = N/L^ and thermodynamic 
potential Qqo 

a = B/Bo {a E (0, 1)), x = p/igBof^ h{x, a) = ^lqol{2g^Bl). (22) 

In the absence of a chromomagnetic field the total thermodynamic potential of fermions 
is equal to 

^qo\B=o = '^{gBofKx,Q). (23) 
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When a chromomagnetic field is generated in the system, B = aB^ ^ 0, the dimensionless 
energy of quarks becomes equal to h{x, a), and the energy of the field itself should also be 
added. Then the total thermodynamic potential becomes: 

52 



^tot — ^qO + — -Sq 



(24) 



Generation of the chromomagnetic field changes the thermodynamic potential of the system 
by the amount 

^ + 2/ KtITX^^ «) - o\s/2 ^ 0) • (25) 



In the model we have chosen, the chromomagnetic field B is a variable parameter and is 
determined by the quantity a, while Bq and p are assumed to be prescribed quantities. 
In this case, a will "adjust" itself in the interval (0,1) with the other parameters being 
fixed such that the total energy of the system takes a minimum value, fitot ^min- Upon 
optimization with respect to a, the dimensionless energy gain due to color ferromagnetism 
becomes a function only of the dimensionless fermion density x: 



X)Ferr = = ^^^^ 



«2 



+ 2g'{h{x,a) - h{x,0)) 



(26) 



The possible generation of a chromomagnetic field and a corresponding ferromagnetic phase 
is in competion with the possible creation of a color superconducting phase in the absence of 
a chromomagnetic field. The system chooses the phase with the lowest energy. The energy 
gain due to creation of a diquark condensate A 7^ can be obtained from equation (fTHl ) 
and is equal to Ssl 

Aficsc = (27) 

Now, to decide which of the phase is preferable, one should compare the energy gains due 
to production of a color ferromagnetic state Afimin and the superconducting state Af2csc- 
For convenience, we introduce dimensionless quantities 

Aeperr = = Acperr ( g^2 (^''5^)3/2 ) , 

— A^CSC _ { A \2,,2 ^ 

Aecsc - —Bf- - -[^) ) 

and compare them with the help of numerical calculations. 

First, the dimensionless energy gain Ah{x, a) = h{x, a) — h{x, 0) was calculated as a 
function of the dimensionless fermion density x at fixed a. It was demonstrated that this 
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difference is always negative, and this indicates that the fermion energy is lowered when 
the field gets nonvanishing values. Moreover the fermion energy oscillates with growing 
density. The reason for this is that Landau levels are filled one by one with growing number 
of fermions. One may expect that at high enough temperature these oscillations should 
vanish. 

Moreover, Ah{x, a) as a function of a dimensionless chromomagnetic field a at fixed 
fermion density was also calculated. It was demonstrated that the stronger the field the 
lower is the energy of fermions until all fermions turn out to be on the lowest Landau level. 
When the fermion density is high enough, this value of a is far above the interval (0, 1). 

In Fig. 1, the total variation of the dimensionless energy of the whole system as a function 
of the dimensionless chromomagnetic field a is depicted. At zero fermion density x = 0, the 
appearance of a nonzero field is not favourable and the dependence is parabolic. At finite 
fermion densities, the growth of the chromomagnetic field is favourable at first, but when 
all the fermions are on the lowest Landau level, the energy begins to grow as a^. For higher 
fermion density, this saturation occurs when a > 1, i.e., the system prefers the favourable 
value of a close to 1 (due to oscillations, not necessarily a = 1 ). 

In Fig. 2, the dimensionless energy gains in color ferromagnetic (with an optimal value of 
a is chosen) and superconducting phases are depicted as functions of the chemical potential 
fjL. The choice of the phase by the system is determined by the condition that the energy 
should take its minimum value in this phase. 

The consideration of the phase transition in the system leads to the conclusion that 
with large enough values of /x, the superconducting phase is preferable. With lowering fi, 
the ferromagnetic phase can become more preferable, and a phase transition takes place. 
However, one may see in the Figure that the energy gain corresponding to the ferromagnetic 
phase is an oscillating function of and hence the picture of phase transitions is somewhat 
more complicated. With further decrease of /i, the superconducting phase may return, then 
again the ferromagnetic phase occurs, and this may repeat several times. 

IV. PHASE TRANSITION AT FINITE TEMPERATURE 

To describe the system at finite temperature, we should use the formulas for the termo- 
dynamic potential of fermions ^q{a, ii,T, Bq), presented in Section 2. Restricting our con- 
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Figure 1: Total change of the dimensionless energy of the whole system as a function of the 
dimensionless chromomagnetic field a. 
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Figure 2: Dimensionless energy gain in a ferromagnetic phase (after a choice of an optimal value 
of a has been done) and in a superconducting phase as a function of the chemical potential [i. 
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sideration only to the main contribution of the "classical " chromomagnetic field -B^/(87r), 
and neglecting small quantum fiuctuations of the gluon field around it, we consider the con- 
tribution of fermions, defined by (|T7jl . (fT^ . where we have to put A = 0, and hence use the 
formula 



Qg = -T^ln 



l + exp( — - — , 



(29) 

fc,: 

It should be mentioned that unlike the case of zero temperature, at finite temperature the 
chemical potential fi becomes an independent variable and is no more equal to the Fermi 
energy. 

As it was done in the case of zero temperature, we employ dimensionless variables and 
quantities, and moreover define a dimensionless temperature 

T 



(30) 



VgBo 

Since \/gBo ~ 10^ MeV, we have for the temperature T ~ 10^ eV ~ lO^^K. 

Now in the case of finite temperature, we also use a numerical calculation method, which 
gives only approximate results, as we have to deal with an infinite series over the fermion 
energy levels. The results are shown in Figs. 3 and 4. 

In Fig. 3, the dimensionless ferromagnetic gain of the total thermodynamic potential, 
minimized with respect to the dimensionless chromomagnetic field a = B / Bq, is depicted: 

[a^ 1 1 
Au;(/i,r) = min„ — + — ^ /x, T, 5o) - (^^(O, /i, T, Eq)) . (31) 

oTT 

Our calculations demonstrate that, with growing temperature, oscillations in the plot of u 
as a function of fi become less evident and, as it is seen in the Figure, in the limit of high 
temperature, they practically disappear. Another interesting observation is that even at 
zero chemical potential the appearance of a finite chromomagnetic field is preferable. This 
is due to the fact that at nonvanishing temperature, fermions exist with finite density even 
at = 0, and the appearance of a chromomagnetic field leads to a finite energy gain. 

Figure 4 shows that with growing temperature and at fixed fi, the thermodynamic po- 
tential gain Au tends to a constant value. It is seen that the higher the chemical potential, 
the lower temperature r(yu) is needed for stabilization of Auj to take place. At low values 
of the chemical potential, Auj decreases with growing temperature and it is stabilized at 
high enough temperature. These facts can be explained in the following way. The discrete 
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Figure 3: Dimensionless thermodynamic potential gain Aw as a function of /x at zero and finite 
(sufficiently high) temperatures. The quantity r is given by equation (30). 
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character of the Landau levels is responsible for the nonmonotonic behavior of Auj{fi). This 
situation does not change while the temperature remains considerably lower than the dis- 
tance between levels in the vicinity of the Fermi surface. It is clear that with growing fi, 
the energy levels become distributed with greater density and hence, lower temperature is 
needed in order to smooth away oscillations. 

V. CONCLUSIONS 

In the present paper we investigated further the gauge field model with a constant chro- 
momagnetic field, i.e., the ferromagnetic state. We demonstrated that the method, proposed 
in [2^, of finding a stabilized solution for this configuration, is valid only if a physically jus- 
tified condition is fulfilled, i.e., the chromomagnetic field exists inside certain domains with 
finite dimensions. This implies that there is a maximum value of the chromomagnetic field 
inside the domain determined by the finite spatial extentnion of the region occupied by the 
field in the direction of the field. 

The main object of our paper was to study the fermion sector of the model, and in 
particular to consider transitions between a color superconducting phase and a possible 
ferromagnetic phase from the point of view of comparing the energy gain in these phases. 
Our observations can be summarized as follows. At comparatively low densities of the 
baryon matter, quarks are confined, and exist only in colorless combinations. With growing 
chemical potential baryons approach each other to such close distancies that quarks become 
free particles. At sufficiently high density, the color superconducting phase can be formed. 
However, for "intermediate " values of the chemical potential, a ferromagnetic phase may 
emerge. There can be a phase transition between these two phases depending on which 
phase has larger energy gain. At zero or comparatively low temperatures a nontrivial phase 
structure can be formed, and this is due to a nonmonotonic dependence of the energy gain 
of fermions in the chromomagnetic field on their density. With growing chemical potential, 
a ferromagnetic state can prove to be favorable, then it is changed by a superconducting 
state, and then again ferromagnetic, and only for high enough chemical potentials, the 
superconducting state becomes always dominating. 

At sufficiently high temperatures (~ ^/gB^), the phase structure is simplified. The 
fermion energy gain is no more dependent on temperature and becomes a monotonic function 
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Figure 4: Stabilization of the thermodynamic potential gain with growing T at various values of 
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of For low ijl the ferromagnetic phase becomes more favorable, and for large enough the 
superconducting phase, there being only one transition point between them. In the temper- 
ature scale, stabilization occurs at T > T(/x), where T{iJ,) is the minimal temperature needed 
to smooth out oscillations of the thermodynamic potential Q that appear due to the discrete 
character of the fermion spectrum in the chromomagnetic field. Theoretical reasoning and 
computer calculations demonstrate that the function T(/x) should be decreasing. 

Further investigations will help to describe more clearly the mechanism of formation of 
the chromomagnetic domains and to consider the boundary conditions in a more sophis- 
ticated manner. Moreover, in order to produce a more realistic picture of the process of 
the ferromagnetic phase formation, quantum fiuctuatious around the "true vacuum " state 
should be taken into account. Estimates showed that at a chemical potential sufficiently 
high for the ferromagnetic state to be formed, interactions between quarks are in no way 
week and should also be considered for. Certain changes to the picture of the process may 
further be given by the interaction of quarks with the charged boson condensate, and this 
should also be considered. Finally, one should also understand that in the case of the SU (3)c 
gauge group, there can be states with a color ferromagnetism in one U{1) subgroup, and 
the color superconductivity in the other U{1) subgroup of the maximal Abelian subgroup 
U{1) X U{1). These problems are to be considered in further investigations. 

Acknowledgements 

Two of the authors (V.Ch.Zh. and O.V.T.) gratefully acknowledge the hospitahty of 
Prof. M. Mueller-Preussker and his colleagues at the particle theory group of the Humboldt 
University extended to them during their stay there. This work was partially supported by 
the Deutsche Forschungsgemeinschaft under contract DFG 436 RUS 113/477/4. 



[1] V.Ch. Zhukovsky, preprint Mad/TH/91-5 (1991), Univ. of Wisconsin-Madison; D. Ebert, 
V.Ch. Zhukovsky and A.S. Vshivtsev, Int. J. Mod. Phys. A13, 1723 (1998) . 

[2] A.V. Averin, A.V. Borisov and V.Ch. Zhukovsky, Z. Phys. C48, 457 (1990). 

[3] D. Ebert and V. N. Pervushin, Teor. Mat. Fiz. 36, 313 (1978); D. Ebert, Yu.L. KaUnovsky, 
L. Miinchow and M.K. Volkov, Int. J. Mod. Phys. A8, 1295 (1993). 



20 



[4] D. Ebert and H. Reinhardt, Nucl. Phys. B271, 188 (1986); D. Ebert, H. Reinhardt and M.K. 

Volkov, Progr. Part. Nucl. Phys. 33, 1 (1994). 
[5] A.S. Vshivtsev, V.Ch. Zhukovsky, K.G. Klimenko and B.V. Magnitsky, Phys. Part. Nucl. 29, 

523 (1998). 

[6] K.G. Klimenko, Teor. Mat. Fiz. 89, 211 (1991); 90, 3 (1992); Z. Phys. C54, 323 (1992); 
I.V. Krive and S.A. Naftuhn, Phys. Rev. D46, 2737 (1992); A.S. Vshivtsev, K.G. Klimenko 
and B.V. Magnitsky, JETP Lett. 62, 283 (1995); Teor. Mat. Fiz. 106, 319 (1996). 

[7] V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, Phys. Rev. Lett. 73, 3499 (1994). 

[8] V.P. Gusynin, Ukrainian J. Phys. 45, 603 (2000). 

[9] K.G. Khmenko, B.V. Magnitsky and A.S. Vshivtsev, Nuovo Cim. A107 (1994) 439; Theor. 

Math. Phys. 101, 1436 (1994); Phys. Atom. Nucl. 57, 2171 (1994). 
[10] D. Ebert and V.Ch. Zhukovsky, Mod. Phys. Lett. A12, 2567 (1997). 
[11] G.K. Savvidy, Phys. Lett., 71B, 133 (1977). 
[12] S.G.Matinyan, G.K.Savvidy, Nucl. Phys., B143, 539 (1978). 
[13] N.K. Nielsen and P. Olesen, Nucl. Phys. B 144, 376 (1978); B 160 380 (1979). 
[14] H.B.Nielsen and M. Ninomiya, Nucl. Phys. B 156, 1 (1979); ibid. B163, 57 (1980); ibid. 

B169 309 (1980); J. Ambjorn and P. Olesen, Nucl.Phys. B170 60, 265 (1980). 
[15] A. Cabo and A.E. Shabad, Trudy FIAN SSSR 111 (1979); A.Cabo, S.Penaranda and 

R.Martinez, Mod. Phys. Lett. AlO 2413, (1995) 
[16] V. Skalozub, M. Bordag, Nucl. Phys., B576, 430 (2000). 
[17] V.V. Skalozub, A.V. Strelchenko, Eur. Phys. J., C33, 105 (2004). 
[18] A.O. Starinets, A.S. Vshivtsev, V.Ch. Zhukovsky, Phys. Lett., 322, 403 (1994) 
[19] S. C. Zhang, H. Hanson and S. Kivelson, Phys. Rev. Lett., 62, 82 (1989) 
[20] A. Iwazaki, O. Morimatsu, Phys. Lett., B571, 61 (2003); |nucI-th/0304005' 



[21] A. Iwazaki, O. Morimatsu, T. Nishikawa, T. Ohtani, Phys. Lett., B579, 347 (2004); 

| hep-ph/ 0309066; Phys. Rev., D47, 034014 (2005); plep^ph/040420'l| 
[22] G. W. SemenofF, Phys. Rev. Lett. 61, 516 (1988). 

[23] M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B422 247 (1998); Nucl. Phys. B537, 
443 (1999) 

[24] K. Langfeld and M. Rho, Nucl. Phys. A660, 475 (1999); |hep^ph/9811227[ 



[25] J. Berges and K. Rajagopal, Nucl. Phys. B538, 215 (1999); pIep^h/98 042331 



21 



[26] T.M. Schwarz, S.P. Klevansky and G. Papp, Phys. Rev., C60, 055205 (1999); |micrt h/9903048| 

[27] M. Alford, Ann. Rev. Nucl. Part. Sci. 51, 131 (2001); |hep^h /0102047| 

[28] M. Alford, Progr. Theor. Phys. SuppL, 153, 1 (2004); [^1101^^0312007, 

[29] LA. Shovkovy, |nucl-th/04l609T] 

[30] D. Ebert, K.G. Klimenko and H. Toki, Phys. Rev. D64, 014038 (2001); pIep^ph/0011273 



[31] D.G. Yakovlev, A.D. Kaminker, P. Haensel, O.Y. Gnedin, Astron.&Astrophys., L24, 389 
(2002). 

[32] F.M. Walter, J.M. Lattimer, Astrophys. J. L145, 576 (2002) 

[33] A. A. Sokolov, I.M.Ternov, V.Ch. Zhukovsky, and A.V. Borisov, Gauge Fields (in Russian), 
Moscow, 1986. 

[34] D. Ebert, V. V. Khudyakov, V. Ch. Zhukovsky, K.G. Klimenko, Phys. Rev., D65, 054024 
(2002); |hep-ph/0106"TlO, V.Ch. Zhukovsky, K.G. Khmenko, V.V. Khudyakov and D.Ebert, 
JETP Lett. 74, 523 (2001). 

[35] K. Rajagopal, F. Wilczek, |hep-ph/0011333l 

[36] Strictly speaking, our following calculations refer to some given macroscopic domain. The 
obtained results turn out to depend on color and rotational (Lorentz) invariant quantities 
only, and are independent of the concrete domain. 

[37] We consider 7— matrices in the 4— dimensional Euclidean space with the metric tensor 
9fiu = diag(— 1, — 1, — 1, — 1), and the relation between the Euclidean and Minkowski time 
^^E) ~ ^^(M)- ^(E) ~ ^')'(A'/)' ^{E) ~ ^(M)- what follows we denote the Euclidean Dirac 
matrices as 7^, suppressing the subscript (E). They have the following basic properties 
7^ = — 7;xj {7/x>7iy} = — 25^j^. The charge conjugation operation for Dirac spinors is defined as 
tpc{x) = C (^(x)"*")* with C^j^C^^ = — 7/x- We choose the standard representation for the 
Dirac matrices (see (24 j). The 75 has the following properties: {7'^,75} = 0, 75^ = 75 = 75- 
Hence, one finds for the charge-conjugation matrix: C = 7*^7^, C"*" = = C* = —C. 

[38] The vanishing of the pion mean-field < vr >= ttq = is here related to the assumed parity 
conservation of the ground state. In what follows we consider only the diquark condensate, 
neglecting the possible existence of a quark condensate < a >= (Tq = 0. 



